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SEM borstI and florian simatos^ 



Abstract. We consider a stochastic network with mobile users in a heavy- traiBc regime. 
We derive the scaling limit of the multi-dimensional queue length process and prove a 
form of spatial state space collapse. The proof exploits a recent result by Lambert and 
Simatos 8 which provides a general principle to establish scaling limits of regenerative 
processes based on the convergence of their excursions. We also prove weak conver- 
gence of the sequences of stationary joint queue length distributions and stationary 
sojourn times. 



We consider a stochastic network with mobile users, originally introduced in Borst et 
al. [2] as a model for a wireless data communication network. Fluid limits of this model 
were studied in Simatos and Tibi |15| and in this paper we examine the heavy- traffic 
characteristics. 

In this model, users arrive at each of the nodes according to independent Poisson 
processes and then move independently of one another while still in service. The tra- 
jectories of the users are Markovian and governed by an irreducible generator matrix 
with stationary distribution n (see Section|2]for notation and definitions) . At each of the 
nodes, users share the total capacity of the node according to the Processor-Sharing dis- 
cipline. This assumption affects the sojourn time distribution but not the distribution 
of the number of users at the various nodes since we will restrict ourselves to the case 
of exponential service requirements. The fundamental difference between this model 
and Jackson networks is that in this model, users move independently of the service 
received: transitions from one node to another are governed by the users themselves 
rather than by completion of service. 
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Coexistence of two time scales and tiie tiomogenization property. The above-mentioned 
feature has an important consequence: although the rate of arrivals to and departures 
from the network is bounded, the rate of movements of users within the network grows 
linearly with the number of users. The model thus shares fundamental characteris- 
tics with two classical queueing models, namely the Ml Mil and the MIMloo queues, 
which are brought about in heavy-traffic conditions. As a result, it inherits the typical 
time scales of both these models: a fast time scale which governs the internal move- 
ments of users between nodes {the MIMloo dynamics) and a slow time scale that gov- 
erns arrivals and departures (the Ml MIX dynamics). In particular, the total number of 
users in the network evolves slowly compared to the speed at which users spread in the 
network, provided the network is highly loaded. 

One of the key technical challenges is to control the M/M/cxD-like dynamics of the 
internal movements of users on the slow time scale of the Ml Ml 1 queue. The main 
idea is that it takes a constant time for any user to be arbitrarily close to the stationary 
distribution n while on the other hand, in a finite time window the number of users 
can only have evolved by a bounded amount because of the Ml Ml 1 dynamics. Thus 
starting from a large initial state, on times of order one the total number of users will 
have essentially stayed the same, while each user will be close to the stationary distri- 
bution n. The law of large numbers therefore suggests that starting from a large initial 
state, after a constant time, users should be spread across the various nodes according 
to n, i.e., approximately a fraction Tif^oi the users should be at node A;. We call this prop- 
erty the homogenization property. To derive the heavy- traffic limit, we need to show that 
the system stays homogenized not only at a given fixed time, but as long as there are a 
large number of users in the network. Technically, this entails control over some hitting 
times, which we achieve via martingale and coupling arguments. 

Convergence of the full process via the convergence of excursions. The homogeniza- 
tion property provides a picture of what happens when there are a large number of users 
in the network: the users are spread across the various nodes according to the station- 
ary distribution n, and in particular it is unlikely for any of the nodes to be empty. Thus, 
the full aggregate service rate is likely to be used, and the total number of users evolves 
as in a single Ml MIX queue with the combined service rate of all nodes. This prop- 
erty provides a useful handle on the processes of interest far away from zero. Imagine 
for instance that the network starts empty: it will eventually become highly loaded, at 
which point the homogenization property kicks in and holds until the network becomes 
empty (or close to) again. In other words, the homogenization property should give us 
control over excursions that reach a certain height and it is therefore natural to expect 
the entire process to converge as well. This line of argument has been used in Lambert 
et al. O to analyze the scaling limit of the Processor-Sharing queue length process and 
has later been generalized in Lambert and Simatos |8|. The proofs of the scaling limit 
results in the present paper leverage the general principle established in Lambert and 
Simatos |8|. 

The above arguments lead to the result that the joint queue length process asymptot- 
ically concentrates on a line whose angle corresponds to the stationary distribution n, 
thus exhibiting a form of state space collapse. The total number of users, after scal- 
ing, behaves asymptotically as in a single Ml Ml 1 queue, and thus evolves as a reflected 
Brownian motion, with the stationary distribution converging to an exponential distri- 
bution. These characteristics are strongly reminiscent of the heavy-traffic behavior of 
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the joint queue length process in various queueing networks, see for instance Bram- 
son [31, Reiman [l3l, Stolyar [16], Verloop at al. [17] and Williams [19]. However, to the 
best of our knowledge, this is the first result which shows that mobility of users, rather 
than scheduling, routing or load balancing, can act as a mechanism producing state 
space collapse. 

Organization of the paper. Section |2] sets up notation used throughout the paper; it 
also presents two key couplings, one with an Ml MIX queue which provides a lower 
bound and one with a closed system that is easier to handle with regard to homogeniza- 
tion. This closed system is analyzed in Section |3| where two bounds are derived: one 
concerning the time needed for the closed system to get homogenized, and the other 
concerning the time that the system stays homogenized, starting from a homogenized 
state. These estimates are used in Section H] to derive corresponding bounds for the 
open system. These bounds allow us to prove that the system is null-recurrent in the 
critical case, a case that had not been treated earlier. 

The last three sections then deal with the heavy-traffic regime: SectionjSjproves, us- 
ing the above excursion arguments, that the sequence of processes converges weakly 
towards a multi-dimensional reflected Brownian motion; Section[6]investigates the as- 
ymptotic behavior of the stationary distributions and Section [7] examines the asymp- 
totic behavior of the sojourn times. 

Relation witli previous work. Some of the results of the present paper can partially 
be found in Simatos and Tibi 1 15 1 and we wish to explain the new contributions. The 
results of Section [3] somewhat strengthen the derivations in flSl : we establish tighter 
bounds and remove a technical condition on the generator matrix governing the mo- 
bility of users. One of the main results of Section ID namely Proposition 14.21 can also 
be found in fTsl under an additional technical condition on the generator matrix. In 
the present paper we provide an alternative proof of this result which we believe to be 
potentially useful in a more general setting. All the other results are completely new. 
In order to have a self-contained paper, and also because we could significantly sim- 
plify some tedious technical detaUs of 1 15 1, we present complete (and simpler) proofs 
for results which were already partially known. 

2. Notation and two useful couplings 

In this paper we deal with multi- dimensional processes, typically taking values in 
IR^ for some K>2, but we also need to consider real-valued processes. It is therefore 
convenient to abuse notation and use the common notation || ■ || to denote the L\ norm 

on every /-dimensional space K'. Thus for each / > 1 and yeW we write ||y|| = |yi|H v 

\yi\. In the sequel, &^{u] for w > denotes a Poisson random variable with parameter u. 
It satisfies the following large -deviation type inequality: 

(1) P{&>{u)>v)<exp{-uh{vlu)),v>u, 

where from now on h{x] - xlogx-i- 1- x. 

2.1. User mobility. In the rest of the paper we fix some integer K>2 and we consider a 
network of K nodes. Let ^ be a Markov process with state space {l,...,K] and generator 
matrix Q - {qte, I < k,£ < K), and let y > be the trace of -Q. We assume that Q is 
irreducible, denote by n its stationary distribution, and define tt = mini<j;<7f :/rj; > and 
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n - maxi<j;<^7rjt < 1. We will need to measure distance to n and introduce the function 
Q-M^ ^ [0, oo) defined by 



with the convention p(y) = if |1 y || - 0. For k- 1, . . . , TC, let be the law of i, started at k. 
For f > 0, let A(r) - maxi<fc,/<jf |Pfc(^(r) - e)-ne\, so that A(f) ^ as t^+oo. For £> 0, 
we define T(e) = sup{r > : A(f) > e}. 

2.2. Sequence of networks. For each n> \, consider (A„,it.l ^ k < K) e [0,oo)^ and 
<k<K)e [0,oo)^, and define A„ = A„,i + --- + Xn.K and /x„ = + --- + pin,K- 
Let Xn - [Xn,k,^ < k < K) he the following cadlag, f\l^-valued stochastic process: for 
f > and k e {\, . . . ,K], x„!c[t) is the number of users at node k at time t in the network 
subject to the following dynamics: 

• users arrive at node k according to a Poisson process with intensity A„ ,fc, and 
arrival streams are independent; 

• users have i.i.d. service requirements, exponentially distributed with parameter 
one; 

• node k serves users according to the Processor-Sharing service discipline and 
has capacity ^in,k\ 

• while still in service, users move independently from everything else according 
to a Markov process with generator matrix Q. 

The memoryless property of the exponential distribution implies that x„ is a Markov 
process with generator Q„ given by 

K K 

o„ (/) (y]^Y.^n,k (/(y + efc) - /(y)) + ^ ii„,k [fiy - sk) - /(y)) Myk>Q] 

k=l k=l 

+ E qkeyk[f{y-ek + ee)- f{y]] 

l<k,e<K 

for any function / : K and any y - [yk, I < k < K) e N^, and where Ck is the kth unit 
vector of M^. Thus x„ can be seen as a system of particles, where particles are added and 
removed (when possible) according to independent Poisson processes attached to each 
node, and while alive move independently according to the same Markovian dynamics. 
Note also that because of the memoryless property, the Processor-Sharing assumption 
has no impact on the law of x„; nonetheless this assumption will impact sojourn time 
distributions studied in Section!?] 

For y E N^, let be the law of x„ started at y (from a network perspective, users 
start with i.i.d. exponential service requirements), and denote by r,j the J^T- dimensional 
process r„ = x„/||x„|| with the convention rnit) - n when ||Xn(f)|l - 0, and where from 
now on if fa E IR and y e then by e denotes the vector {byk, l<k<K). 

We denote by a^ k the arrival process at the fcth node and by d„^k the potential depar- 
ture process from the fcth node, so that a^^k is a Poisson process with intensity A„ and 
d„^k is a Poisson process with intensity idn,k such that the 2K processes {an,k>dfi,k, 1 ^ 
k< K) are independent. Moreover, by definition of Xn it holds that 

K K r- 

(2) l|x„(f)|| = |U„(0)||+Ean,fc(f)-L / 

{u-)>0}dn kidu), t>0. 

k=l k=lJl0.t] 

Let an - an,\ H v an,K and d„ - dn,i H H dn,K, so that a„ and d„ are independent 

Poisson processes with intensities A„ and respectively. Define p„ - An/ fin- it has 
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been proved in |5l[T5] that Xn is positive-recurrent when pn < 1 and transient when p„ > 
1; in Proposition 14.31 we will complete the picture and prove that x„ is null-recurrent 
when pn - I. When Pr, < 1 denote by v„ the stationary distribution of x„. We define 
K - sup„(A„ -I- fin) and assume throughout the paper that k is finite. 

2.3. Functional operators. Fix some / > 1 and £ > 0. Let Di be the space of IR' -valued 
cadlag functions. For / - {fjc, I < k < i) e Di, we define the following operators: 

rk/,£) = inf{f > 0: ||/(f)|| < £}, Thf,E) = inf{f > : ||/(f)|| > e] 

as well as 

ro(/) = inf{?>0:||/(r)||=0} and fo(/) = min Toift). 

l<fc<! 

For t> 0, let a,6t : Di Di be the stopping and shift operators, defined by a{f) (s) - 
/(sA Toif)) and dtif){s) - f{t+ s), respectively, for / e Di and s > 0. Define also the 
map e| : D,- Di as follows: 

4(/) = (f7°0ri(/,.))(/)- 

In words, el (/) is the process / shifted at the first time r' (/, e) when ||/|| reaches level e 
and stopped at the first time it reaches afterwards. Finally, let (/) be the left endpoint 
of the first excursion of ||/|| to reach level e: 

g,(/) = sup{r<r'(/,£):||/(flll=o}. 

Note that similarly as ||-||, we use the same notation to refer to operators defined on 
functions taking values in IR' for any / > 1; for instance, we have (/,£) - r' (||/||,£). 

2.4. Coupling with an Ml Mil queue. The process Xn is naturally coupled with the 
queue length process of an M/M/1 queue with arrival rate A„ and service rate p„. First, 
for any f eD\, define the function /, called the function / reflected above its past infi- 
mum, by 

/(f) = /(f) -min inf /(s),0 . 

— \0<s<f j 

For n>\, let in the sequel - ||x„(0)|| -i- a„ - dn and - Ln process £n re- 

flected above its past infimum. Since fl„ and d„ are independent Poisson processes with 
respective intensity A„ and [in, (n is a continuous-time random walk and (n is equal in 
distribution to the queue length process of an Ml Ml 1 queue with arrival rate A„ and 
departure rate p„. 

Intuitively, this coupling does the following: the potential total departure process 
from both x„ and is given by d„. When d„ jt rings, there is no departure if x„ ^ - 
while there may be other users elsewhere. The process (yi ignores how users are spread 
in the network: if there are users in the network and one of the d„ jt rings, then one of 
the users leaves. Thus || x„ || and (n coincide as long as there is no empty node, and at all 
times there are more departures from (yi than from x„. Formally, we have the following 
result; recall that if p„ < 1, then x„ is positive-recurrent with stationary distribution v„. 

Lemma2.1. Foranyye N^, the two following properties holds P^''^^'^ost surely: 
. ||x„(f)|| = ^„(f) /or a// f< fo(x„); 
• \\Xn[t)\\> £nit) forall t>0. 

In particular, if pn < 1, then ||x„(0)|| underP^' is stochastically lower bounded by a geo- 
metric random variable with parameter pn- 
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Proof. We prove the two first properties by induction. By construction, it liolds that 
IUr((0)|l - ^n(O) and when a„ rings, both processes ||x„|| and £n increase by 1, which 
preserves the difference ||x„|| - f „. Internal movements in x„ also preserve this differ- 
ence, hence one only needs to inspect what happens when one of the rf„ ^ rings. 

So to prove the first property, consider t<To (x,,) and l< k< K such that rfn.fc (i t]) - 1 
and assume that ||x„(f-)|| - £n{t-): we must show that ||x„(f)|| - (n{t). Since t < 
ro(x„), by definition of ToCxn) we have Xn,kit-) > and so ||x„(t)|| = ||x„(f)|| - 1. On 
the other hand, we have £n{t-) - II Xf, ( f-) || > by induction hypotheses and so by con- 
struction, ( nit) - a nW - 1, which proves the desired property. 

Let us now prove the second property, so consider f > and 1 < A; < such that 
dn,ki{t]) - 1 and assume that ||Xrt(r-)|| > £„{t-]: we must show that ||x„(f)|| > £n[t]. If 
£n[t-) - then ^„(r) - and ||x„(f)|| > £„[t). Else, £„[t) - £„[t-) - 1 and since ||x„(f)|| 
decreases by at most 1 (it does when x„ j;(t-) > 0, otherwise it stays constant) we also 
have ||Xn(t)|| > ^„(f) in this case. 

As for the second assertion of the lemma, for any q >0 we have by stationarity and 
using the second property 

Pi" {\\x„m\\ >q]^ P^" (l|x„(f)|| > ^7) > Pi" [£n[t) > q] — (pj"? 

'+00 

since £„ ( f) converges in distribution as t goes to infinity to a geometric random variable 
with parameter p„. This proves the result. □ 

We will often use the previous lemma in combination with the following lower bound 
on fo(x„): 

(3) Thrn-n,6)<fQ[Xn), S<n. 

Indeed, since by definition we have \\rn{t) - 7r|| > |r„jt(f) ~Ttk\ we see that Xn,k{t) - 
implies that ||r„(f) - n\\ > n. In particular, this implies that ||x„(?)|| = £nit] for all 
t< [rn-Ti,8) and S <n. 

2.5. Coupling with a closed system. Let xj, be the process built on the same proba- 
bility space as x„, sharing the same stochastic primitives as x„ but ignoring arrivals 
and departures. More precisely, if ((n,i>^ - i - Ilyll) are the ||y|| independent (but not 
identically distributed, due to the initial conditions) trajectories of the ||y|| initial users 
under P^, then we define 

Ilyll 

= £ W=fcl' t>0,l<k<K. 
1=1 

Note that x^ is a Markov process with generator O' defined similarly as Qn but with all 
A„ fc's and /x„,fc's equal to 0. In particular, the law of x'^ does not depend on n. Under P^, 
by construction it holds that of the {^n,i,i ^ i ^ Ilyll) are i.i.d. with common distribu- 
tion ^ under P^. 

We define r'„ - x'„l\\x'„\\ with the usual convention r'„{t) -n when ||x5j(?)|| - 0. Note 
that because xj, is a closed system, we have || x'^ {t)\\ - \\ x'^ (0) || . The following inequalities 
are intuitively clear, see for instance 1 15 1 for a formal proof: for any t>Q and k- \,...,K, 

(4) - dnit) < Xn,k{t) - x'^ ,^{t) < Unit) and -dn{t)<\\xn{t)\\-\\x'n{t)\\<an{t). 
This has the following useful consequence. 
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Lemma 2.2. For any ye N with \\y\\ > 0, we have 

vl , 2K{an[t] + d„(t)]\ 

K (Vf > : ||r„(f) - r'jm < "^^^^ = 1. 

In particular, for any y e N^, any 5 >0 and any t>0, we have 

(5) py„{\\rn{t)-r'^{t)\\>S]<P{&>{Kt)>S\\y\\/(.2K)]. 

Proof. Fixk^l,...,K and yeN^ with ||y|| > 0: then under P^, it holds that 



\r'k^t)-rn,M\ 



\\Xn{t)\\ 



x'^l,{t)\\Xn{t)\\-Xn,k{t)\\ 



llyll IU„(f)li 



ix'^j^it) - X„_k ( ?)) II Xn ( f) II + Xn,ki f) (II X„ ( f) || - || y ||) 



llyllllx„(f)|| 

K^t^t)-x„,k{t)\ ^ |||x„(f)||-||y||| 



Together with (4) this gives the first result, which implies |5) since an{t) + dn( t) for 
any f > is a Poisson random variable with parameter (A„ + t and A„ + < jc by 
definition of jc. □ 

3. Control of the closed system 

In this section we are interested in the closed system xj,. Since its law does not de- 
pend on n, in order to simplify the notation we remove temporarily aU subscripts n and 
write x' and r' instead of Pj^, and r^. We will denote x'{t) = l{f,.(o=fc!, with 

{^i) independent Markov processes with generator matrix Q. 



3.1. Homogenization at a fixed deterministic time. We first show that starting from 
any initial state, the system becomes close to homogenization in a constant time. Note 
that the following bound is consistent with the central limit theorem, which suggests 
that \\r'{t)-n\\ should be of order ( || x' (0) || ) " ^ for large f and || x' (0) || . The following 
result improves on Simatos and Tibi (TS] Proposition 5.2] ; here we use Chernoff 's instead 
of Chebyshev's inequality. 

Lemma 3.1. There exists eq > 0, depending only on n and K, such that for any Q<e<eq 
and anyy e N^, 

/ f2|, 

py (II r'[T{El[2K))) - :fr|| > £) < 2ii:exp -- 



\ 4^2 

Proof In the rest of the proof, fix £ > 0, y £ N^, and write t-T[El(2K]) and£' - eI[2K). 
Since r'(f) -n under P°, the bound holds when ||y|| - and we consider ||y|| > 1. Stan- 
dard manipulations yield 



P>'{||r'(f)-:?r||>£]<^: max ff 

i<fc<jf 



\r'^{t)-nk\^ElK]. 



^ llyll 



Until the end of the proof fix some 1 < A; < ^C: we have x^(f) - l{f,(o=fc) 
and so 



=ic\ under P^ 



\£Hr'^{t]]-nk\^ 



I llyll 

y py{£,i{t)^k)-ni 

llyll a 



1 llyll 
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Thus by definition of A and t, we have |E^(r^(f)) - Jijd < e' since t - T(e'). Conse- 
quendy, die triangular inequality gives \ r'j^{t) - n k\ < \ r'j^{t) - Ey {r'j^{t)) \ + e' and so 

py{\r'i^{t)-nk\ > £/K) < P^(|r;(f) - E^(r;(f))| > e') = FJ'(|x;t(f) - E^(x;(f))| > £'||y||) . 

Let us now define £o- For p e [0, 1], define f{p) - p{l - p). Since /(I - tF) > 0, there 
exists £g < 1 - which only depends on'jf, such that / (I - tt- 5) > 5 for all 5 < £g. We fix 
such an e'q and consider eq - 2K£q, which only depends on n and K. In the sequel we 
assume that £ < £o, or equivalently, £' < £g. 

Recall Chernoff's inequality: if (7,-, 1 < i < I) are independent random variables with 
I F/l < 1 and E( Yj) - for each !</</, then for any < 77 < 



i=l 



: rjb 



:2e"''''* with = E(Yf) + --- + E(F/). 



Denote pi - lP^(fi(f) - k): we wish to apply Chernoff's inequality to the random 
variables (7;,! < i < ||y||) with 7,- = l{f,(t)=fc} - Pi, for which fo^ = /(pi) + ■■■ + /(pnyn). In 
order to ease the notation, we suppress the dependencies of pi, Yi and b on k, t and y, 
which have been fixed once and for all earlier. Define now 77 = £'||y||/fo, so that 



llyll 

i=l 



: rjb 



Assume for a moment that ri<b: then we could apply Chernoff's inequality and get 



lljll 
i=i 



: rjb 



< 2e-r'' = 2 



exp 



(e'llyll)' 



< 2exp 



4K2 



using b 



II /4. This would prove the result, and so it remains only to prove that T]<b, 



or equivalently, b'^ > e' 



'2 > £' llyll. By definition we have b^ - f{p\) + — v /(Pllyii)- Let 1 < j < 
llyll: since f > t{e'), we have \pi-nk\ <e' and in particular 

n_- e' <nic-E'< Pi <njc + £' < tt + £'. 

Since f{p) - /(I - p], f is increasing on [0, 1/2] and decreasing on [1/2, 1] and n < 
1 - TT, the previous inequalities imply that 

/(pi) > min (/(TT -£'), /(7f + £')) = min (/(TT -£'), /(I - 7f -£')) = /(I - TT -£')> £' 

by choice of £q and since f' < £q . Thus fo^ > £' || y || which concludes the proof. □ 

3.2. Deviation time from the equilibrium. We now study the time needed for the pro- 
cess r' to leave a neighborhood of n. Note that x' can be seen as a multi- dimensional 
Ehrenfest urn, where sharp results on hitting times in the two-dimensional case K - 2 
have been established in Feuillet and Robert |4|. 

Estimates on hitting times will follow from the optional sampling theorem applied to 
the martingale of Proposition l3.4l This martingale was first constructed in Simatos and 
Tibi 1 15 1 for the open system under an additional diagonalizability assumption on Q. 
Since the martingale construction is quite complicated for the open system, we adopt 
here a different approach: we only construct the martingale for the closed system and 
then use coupling arguments to transfer results on the closed system to the open one. 

Compared to Simatos and Tibi 1 15 1, the new contribution of the following construc- 
tion consists of Lemma \32[ which makes it possible to drop the diagonalizability as- 
sumption. Nonetheless, the construction of the martingale for the open system is in- 
tricate while it becomes quite elementary for the closed one. For this reason, we have 
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chosen not to refer to US], but rather to provide a self-contained proof (except for the 
proof of Lemma l373l which can be repeated almost verbatim) . 

3.2.1. Additional notation. The following notation holds throughout the rest of this sec- 
tion. Let 5^ = {M £ [0, : || u\\ < 1} c K^"! and = {u £ [0, 1]^ : || u\\ = 1} c be the 
^C- dimensional simplex. Let L : 5^ — ► S^k be the function that completes u e into a 
probability distribution, i.e., (Lu)fc -Ukifl<k<K-l and (Lu)k = 1 - II u|| for ue ^. 
Note that L is invertible with inverse : being the projection of the K- \ first 
coordinates. Let finally n - dL\a.%{n\,...,nK) be the diagonal matrix with entries (:/rfc) on 
the diagonal. 

Let / be the Jordan normal form corresponding to Q with change of basis matrix w. 
Thus / and w are possibly complex matrices, and we have the following properties, see 
for instance Herstein and Winter [6|. Let \ < i < I) be the / distinct eigenvalues 
of Q with {t] = 0, for 1 < / < / let m, be the algebraic multiplicity of di and let ci)kii) 
for some k{i) e {l,...,K} be any eigenvector of / corresponding to the eigenvalue di, 
i.e., Jci)jc(i) - diCJkii)- Since m,- is the algebraic multiplicity of di we have in particular 

midi H h mj-idj-i - -y (recall that y > is the trace of -Q). Moreover, we have 

Q = (i)~^Ja) so that e^*? - cj'^e^^oj for any f £ R, and because of the block structure of / 
this gives {e'-^) k[i),j - for j ^ fc(0 and {e'^^] k(i),k{i) - ■ ^he sequel we consider the 
following function F : [0, oo) : 

/-I 

F{u)^l\\{wu)k(i)\'"\ ueU^. 

i=l 

Then F satisfies the following simple property, which is key to generalize the martin- 
gale construction of Simatos and Tibi [15] to the case of non-diagonalizable Q. 

Lemma 3.2. For any u £ IR^ and teU, we haveF[e*'^u) - e~'^*F[u). 

Proof. We have e*'^ - u)~^e*^(i) so that 

i=l (=1 (=1 

The identity {e''^ (t)u)k{i) - e'^' (ww)fc(i) follows from the fact that (e'-')fc(,)_j is equal to 
for i ^ k{i) and to e'-^' for j - k{i]. This proves the result. □ 

3.2.2. Upper bound on hitting times. Proposition l3.5l is the main technical result of this 
section, which will be used in the following section. We omit the proof of the following 
result, for which one can repeat almost verbatim the proof of Lemma A.5 in Simatos and 
Tibi 1151. 

Lemma3.3. The quantity supQ^^^i[c^ f^iP{n~^Lu))'^~^du) is finite. 
Proposition 3.4. For any c> and any y e N^, the process 

McW^e-'rtf f\i^J^Y'^'\F{U-'Lu)Y~' du, t>0, 

is a bounded martingale under . 

Proof. Fix in the rest of the proof some c> and y £ . Then under P^, we have for 
any f > 

Mc(r)<7r-"J'" f [F{n-^Lu)Y~^ du 
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and so supf>QMc(f) is bounded by Lemma 13^ for < c < 1 while for c > 1 one only 
needs to use the fact that uEr^ ^ F{U~^Lu) is bounded. Let s, f > 0: we have 



-cy[t+s) 



r ^ 



exp 



/ K 

U=i 



U'(f) 



with Gfc(M) = log((LM)fc/:/rfc),sothat e*^*'"' = {Yr^Lu)^. For any z £ N^, we have 



exp 



,k=l 















exp 








V 


ll^ll 







M=fc} 



1=1 



exp 



I K 

H Gfc(u)l{f,.(5)=fc} 



since the < i < \\z\\) under are independent. For j e {l,...,K}, zj of the (f/) are 
i.i.d. with distribution ^ under Pj and so 



exp 



' K 

H Gs;(U)l{f,.(s) = fc) 

Vfc=l 













exp 









£ Gkiun 



{f(s)=fcl 



K I K 

n E 

j=iVfc=i 

K I K 



;=i Vfc=i 



Vl K 

n 



Thus 



■'^ fc=i 

We want to make the change of variables iLv)ic/7ik - ie^'^U~^ Lu)^, i.e., Yl~^Lv - 
e^^U'^Lu or Lv = Ue^'^U'^ Lu. To use the inverse of L, we check that Ue^'^Yl'^ Lu e Sfi}^ 
for ueS^: we have {Tle^^Tl~^ Lu] s; > because Li/ has only positive coordinates and the 
matrix ne*'3n~^ has only positive coefficients, and 



X (ne^'3n-iLu)fc = E ^fc E — - = E E ^fcPfc(f(^) = 7) 



E 

jfc=l 



fc=l ;=1 



j=l ^'J fc=l 



= E E(-^")fc=l- 

fc=i 



Hence we can consider v - HgU with /fj - L~^Xle^'^Y\~^ L. Clearly Hs is invertible 
with inverse H'"^ = L~'^Ile~^'^Yr^L = H_s. For any f e K, we have (^) c 5^: indeed, 
HfU for ue ^ has only non-negative coordinates, because of the same arguments as 
above, and moreover 

K-l K-l K-\ K 

fc=l fc=l fc=l fc=l 

Hence //s(^) c ^ and H-A^) = (^) c ^ and so the restriction Hs : ^ 5^ of 
to 6^ is invertible with inverse H-s. This gives 



>f n(^ 



(F(n-iLH_,i;))^-i|Jac„(H_s)|dy 
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where if M is a matrbc |M| stands for its determinant and if M -.W ^ then JaCy(M} 
stands for its Jacobian matrix evaluated at v. 

JaCi,(M] = i^^[v], \<k<b,\<e <a 
\ du( 

Lemma l3^ therefore gives, since Y{~^LH-sV - e'^'^rr^Ly, 

We now show that |JaCy(H_s)| - e^^ , which will complete the proof. The chain rule 
gives 

Jac„(J/_s) = JacL„(L-ine-''?n-i)Jac„(L) = JaCi„(ne-*'?n-i)Jac„(L) 

using for the second equality that is the projection on the first ^-1 coordinates. 
Also, Jac„(L) is the diagonal matrix diag(l, 1, . . . , 1, -1), so its determinant is equal to -1 
and on the other hand, if M is a linear operator then Jac^ (Af) = M for all v and so 

|Jaci„(ne"''?n-i)| = \Yle~'^n-\ = \e~'^\ = e^'. 

The proof is complete. □ 

Proposition 3.5. There exists a family of finite constants (^^§,5 > 0) such that for every 
t> l/y, every d > and every y e with g{y] < nS'^/Q, 

P^(r'(r'-7r,<5) < t^<'^sexp{Klogt-5^\\y\\l8). 

Proof. Fix in the rest of the proof S > 0, f > l/y and y e N^, and denote c - 1/(7?) and 
r = {r' - 71,6): then Markov inequality gives 

py (T<t)< e^y (exp [-cjT]) . 

We now derive an upper bound on this Laplace transform. Consider the martin- 
gale Mc of Proposition l3.4l Since under P'' it is a bounded martingale by Proposition l3.4l 
the optional stopping theorem gives 

E>'(M<:(r)) = EJ'(M,(0)). 

We will provide an upper bound on (Mc (0)) and a lower bound on E^ (Mc (T)) of the 
form j4E^(e~'^^^), thus providing a desired upper bound on E^(e~'^l'^). For u,v e let 
H[u, v) - ui log(ui lvi) + — yuK\o%{uKl vk) be the relative entropy, so that Mc under P^ 
can be rewritten as 

Mc(r) = e-'"^' \ ex:p{\\y\\[H{r\t),n)-H[r'[t),Lu))}[F{n-'^Lu)Y~^ du. 

The relative entropy is positive and satisfies the following upper and lower bounds: 

1 p 1 

(6) -^u-vT <H(u,v)< u,veS'K- 

2 mini; f fc 

The lower bound, called Pinsker's inequality, is well-known, see for instance Pinsker [TT| , 
while the upper bound follows by convexity: 

H{u, i;) = ^ Mfc log 1 -H < ^ < ^ < — , 

fc=i V Vk I fcfi I'k fc=i Vk mmfcVfc 
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using ujc<l. Thus using that the relative entropy is always positive, we obtain using the 
upper bound in (6) 

Ey [Mem) < ell^««<3'/llyll.-) f iFin-^Lu))'-' du < ^igll^llf W'^c"^ = A^eKlogt^yWeiyVn 
iff" 

with A\ - supo<c<i(c^/5o(F(n~^Lu))'^~^(iu) and A2 - Aij^, Ai and therefore A2 being 
finite by Lemma l33l frecall that < c < 1 by assumption). We now derive a lower bound 
onMc(r). The lower bound in (6) gives H{r'{T),n)> ||x(r) - :?r||2/2 and since ||x(r)- 
n\\ > 5 by definition of T we obtain the following lower bound on the integral part of 
Mc[T): 

[ ^\\ymir'm.n)-mr'm.Lu)) (p^n-'mY-'du 

By definition, F is a continuous function, hence it is bounded on the compact set 
^' = n-iL(^)and so 

with A-i^{\ + sup_s^,F)-i. For v ^ S^k let Sfi{v) ^{ue5^K- H{v, u) < S^/A] and cpiv) = 
/^(j^j du be the volume of S^{v): then 

with Ai - infyg^j, (p{v). We have ^4 > since (p is easily seen to be continuous, (p(.v)> 
for any v e and is compact. Therefore, we have proved 

MciT) > e-'rTe\\y\\s'/2^^^-\\yi\s'i4^ ^ ^^^-cjT ^\\y\\S',4 
with A5 - A-^AicL deterministic constant. Combining all the previous bounds, we obtain 

\py{Thr'-n,S)<t]<^exp[Klogt+\\y\\g{y)/n-\\y\\6^/4). 
^ 'As 

Assuming p(y) < nS^/S ends the proof, since A2 and ^5 only depend on 5. □ 

4. Control of the open system 

We extend Lemma [TTI and Proposition 13.51 to the open system, exploiting the cou- 
pling with xjj of Section [231 In the rest of the paper, £0 denotes the constant given by 
Lemma lsTl Recall that k - sup„>j (A„ + is assumed to be finite. 

Lemma 4.1. For any w > 1, any 6 < Ieq and any y e , we have 
(II rniT] - :^|| > 5) < P [Sfi(KT) > 5||y|| /(4ii:)) + 2^:exp 

where T- t{8I(AK)) 
Proof. Since 

K (llrnW - ;r|| > 5) < (II rniT) - r'„[T)\\ > 812) + P^ (|| r^W - ;r|| > 812) 
the result follows directly from Lemmas l2.2l and l3TT] □ 
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Proposition 4.2. There exists a finite constant ci > and for each S > Q, there exists a 
finite constant cz (5) such that ci and cz [S) only depend on K, k and Q (and d for cz (5) ) 
and such that for every n>l, every f > and every <S <l,<peN and y e such that: 

7]<2eo, > max |t (77/(4^0 ||y|| and g{y]<ri, 

whererj - nd^ 132, then 

Pi[Thrn-n,8)<tAT^ (^)) < cz (<5) exp (log t-ci[S*4>- log0)) . 

Proof. In the rest of the proof, fix n, t, 5, (p, rj and y as in the statement of the proposition 
and denote for simplicity t - T(j]li4K)) and u - cprjliSKK). We have 

Lt/uJ-l 

Pi{Tr<tAT^)< Pi{Tr<Tj„iu<Tr<ii+l)u) 
i=0 

where from now on T,- - T' (r„ - 7i,d) and - {x„,(p). Since ||y|| > (p and Q{y) < rj, 
the term corresponding to / = in the above sum is upper bounded by 

¥-i{Tr<T^,Tr<U)< sup wiiTr^U). 

y':\\y'\\>Lip,Q{y')^r) 

Consider now / > 0, and note that / u > t by assumption. Since ||x„(/m- t)\\> <p in 
the event {iu< Tx), the Markov property at time iu-T gives 

Pi[Tr<Tjc,iu<Tr<ii+l)u, II r„ (i M) - ;r|| > 77) < sup p( [\\ rnW - ;r|| > 77) . 

y'.\\y'\\r^(l> 

Similarly, since ||x„(/m)|| > (p in the event {iu < Tx], the Markov property at time iu 
gives 

Pl[Tr<Tx,iu<Tr<U+l)u,\\rn(iu]-7i\\<r])< sup p{(Tr<u). 

y':|l>''||>0,e(y')<»; 

Since the previous upper bounds do not depend on i, summing over < i < [t/ u\ -1 
gives 

/ , , ] 

Pl[Tr<tATx)<[tlu\ sup (||r„(T)-;r||>77)+ sup [T,- < u] 

\y':\\y'\\><P y':||/||>0,e(y')<77 

Let y' e with ||y'|| > (f>: since 77 < 2co and t - TirjUAK]), Lemma l4?T] gives 



P( (l|r„(T) -;r|| > 77) <P[g^[KT) > T](p/{4K)) + 2Kexpl^-^^^ . 
Since cp > SkKt/ti, we have 



P[&>iKT) > ri(l)/{4K)) <P{&>i(l)T]/{8K)) > (pT]l{4K)) < exp 
and we finally get 



y'.\\y'\\^<P 



8K J 



sup P^' (II r„(T) - ;r|| > 77) < exp (" ^^|^) + 2^ exp (-^) ■ 



Consider now y' e N with ||y'|| > (p and p(y') < 77: then 

P^' [Tr<u]< P^n [r' (r„ - r^, 5/2) < j/j + p{ [r' (r^ - ;r, 5/2) < u] . 
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On the one hand, Lemma lZ2l implies that 

(r'(r„ - C6I2) < m) < P( (r'(a„ + rf„,5||y' ||/(4^0) < <¥[&'(ku) > J)\\y'\\im)] 

using that 5 > 77 and that a„ + is an increasing process to get the last inequality. Since 
II y' II > this gives 

[t^ irn - r'n,S/2) < u] < P(3«(7cw) > 84>K4K)) < exp 

plugging in the definition of u and using (1) for the last inequality. On the other hand, 
since u > I/7 and giy') <r]- 7r(5/2)^/8, Proposition l3 .5 l implies that 

'P'n[Thr'n-n,5l2) < u] < <^5/2 exp (riogM- 5^11 y'||/32) 

S '^5/2(77/ (8^:?c))^ exp (^:iog0- 5^0/32) 

using II y' II ><p and the definition of u to obtain the last inequality. Gathering the previ- 
ous upper bounds, we have proved at this point that 

^P^ ( r' (r„ - ;r, 5) < f A rkx„, 0)) < 2^:exp (- ^] + 2 exp (- ^^1^] 

+ '^812 (J?/ (^Kk))^ exp (^:iog0 -8^(pl?,2). 
The result then follows easily from this expression. □ 

Proposition l4.2l has the following consequence, which is interesting in its own right 
and will be used in the proof of Lemma |53] 

Proposition 4.3. Letn> 1; ifpn - 1 then x„ is null-recurrent. 

Proof. Since ^„ is null-recurrent it is enough to show that x„ is recurrent by Lemma [2?T] 
Intuitively, the drift of x„ is always strictly positive due to the fact that there is always 
a positive probability that a potential departure finds an empty node, creating a slack 
between the arrival and departure rates. However, the drift should go to as the size 
of the initial state becomes large. Theorem 3.2 in Lamperti 1 10 1 asserts that if the drift 
vanishes sufficiently fast, then x„ is recurrent. 

Since n is fixed, in order to ease notation, we omit the subscript n, so for instance a 
and d refer to the processes and dn, respectively, and x refers to %„. Let b- 1/4 and 
(p(y) = ||y||i+2b ^ ||y||3/2 fQj. y £ pqj. I > ^ Yet coi = cji-i + (p{x[a)i-i]) with wq = and 
Xi - If we can show that 

(7) ess sup E°{xi+i-Xi\Xi^^'X]^ijJ ^i- l) 

< ^ess inf E° {Hm - Xif I Xi = ^,Xj = (jJ <i-l] + Cr''" 

for some finite constant C and all ^ large enough, where the sup and the inf are taken 
over i > 1 and {^j,0 < j < i - 1), then Theorem 3.2 in Lamperti fTO| will imply that x 
almost surely visits infinitely often some finite set; since it is irreducible this wiU prove 
recurrence. Let - (T{Xj,0 £ 7 £ i) and 'S^,- = a{xiaij),0 < j < /): we have 

E°[xi+i-Xi\^i] = E° [e" iXi+i-Xil'Si) = E° [E^"^'-' {\\x{ip(xmm - llx(O)ll) 
and hence 

E%Xi+i-Xi\^i)^ max (||x((^3,)|| - ||x(0)||) , 

y-\\y\\=Xi 
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writing indifferently q}y or q>{y). Similarly, 



E°{{XM-Xif\^i)> inf Ey[[\\x{(pyn-\\xm\\f). 

y-\\y\\=Xi 



Remember that £ - ||x(0)|| + a - d: since ||x|| > ^ by Lemma \2A\ and £ - £> £ hy 
definition, we get (||x((py)|| - |U(0)||)2 > cpy in the event {£((py] - ||x(0)|| > cpy^} and so 

Ey[i\\x[(py)\\ - iix(0)ii)2) > Ey[[\\x[(py)\\ - \\xmf;h(Py) - umw > ip],'^] > ccpy 

with c - inf (>o P° {£[t]> \/t] which is strictly positive since £{t)l\/l under P° converges 
in distribution as ? ^ +00 to a normal random variable. Thus to show Q it is enough to 
show that there exists some finite constant C such that 

(8) max ( || x{cp^) \\ - || x(0) || ) < 

for all ^ large enough. Let f > and y £ such that \\y\\ integrating (2) over and 
using X-n gives 

Ey[\\x{^>l_)\\-\\x(m)^Y.^^k [XkW ^ o) du. 

k=l Jo 

Letke{l,...,K}: then 
jj' py (Xfc(M) ^0]du< (p^py\^T^ (X, ^W2) < (p^^ + jj'^ py [xfc(u) = 0, r' (x,^/2) > ] dM. 

Since £ < || x|| and ^ is a symmetric random walk, we have 

(r' (x,^/2) < (p,.] < P>' [r^ {£,(12) < = P" (r^ i£,^/2] < (^^ ] = 2P'' (/((^),.) > ^/2) 

using the reflection principle to get the last equality. Moreover, writing for f > as 
£{t) - {a{t) - At) + {At - dit)), using the triangular inequality and the fact that a and d 
are i.i.d. Poisson processes with intensity A, we obtain 

P° {£((pO > ^/2) < P (|i^(A(p^) -A(p^\> ^/4) . 

Extending (1) to upper bound Pi&'iAcp^) < q)(-(/4), it can be proved that there exists 
a finite constant H>0 such that 

P [\&>{A(p() -A(p(\> ^W4) < exp [-m^/q)^] = exp [-H^/fj . 

Further, let 5 > and rj - 7r5^/32 be such that 5 < n_ and 77 < 2eq, and write t = 
T{r]l{AK)]: then 

l^"^' P^ (xfc(M) = 0, r' (X, ^"/2) > (pfj < T + P>' [xfc(u) = 0, r' (x,^/2) > (p,.] 

and for u>t, the Markov property at time t gives 

P-" (xfc(M) = 0, (X, ^/2) > < P-" (II r(T) - ;r|| > 77) 

+ sup P^' [xfc(u - T) = 0, (X, ^/2) > (^^c - rj 
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where the supremum is taken over the set {y' e : |iy'|| > £,l2,Q{y'] < tj}. Since t] < lea 
and T - T{ri/{4K)), Lemma fOI gives 

(II r(T) -7i\\>rj)<P [&>{KT) > 6£,im)) + 2Kexp [-^^^ 



; exp + 2_*:exp 



4Kkt ) [ 16K^ I 

using (D for the last inequality, together with h{v) >. v for v large enough (so this in- 
equality holds for ^ large enough). Since Xk{u) -Q implies || r(M) - :7r|| > n_> 5, we have 
for any t < u<ip^ and any y' £ with || y' II > ^' 12 and p(y') < r] 

P^'[xfc(w-T) = 0, rUx,^/2) >(/)^-Tj <P>''[rUr-7r,(5) < w- t, rkx,^W2) > (/)^ - rj 

<V>y' {T\r-n,8]< cp^ A (x, ^/2)] 

< C2((5) exp (log(^f - ci ((5*^/2 - log(^/2))) 



since all the assumptions of Proposition l4.2l are satisfied (for ^ large enough). Gathering 
all the previous bounds, we obtain for ^ large enough 



(II x((p^] II - II x(0) II) < 4^i(pf exp (-H V^) + + p exp (" 



+ 2Kfiexp [- + C2{S)i^(p^ exp (log^^ - ci ((5*^/2 - log(^ /2))) 

which proves J8) since every term decays exponentially fast, except for the constant 
term Tfx which has to be compared to cp^l^^^^ which goes to infinity. □ 

5. Scaling limit 

Fix from now on two parameters A > and a > 0. In the rest of the paper, we assume 
that the following assumption holds. 

Heavy-traffic assumption. From now on, we assume that Pn^^ for each w > 1 and that 
lim An - A and lim n{l - pn) = cc. 

We believe that the techniques of the paper could be adapted to the case a eR (and 
hence remove the assumption pn < 1). Under the heavy- traffic assumption, we have 
jj-n A. Note that we do not require each A„ j; or j^^.k to converge, but only the corre- 
sponding sum. 

In this section we are interested in the sequence (X„, n > 1) of renormalized pro- 
cesses where X„ for each n > 1 is given by 

Xnin^ t) 
Xn{t)^— -, t>0. 

n 

We define Rnit) - r„[n^t) which satisfies i?„(r) - Xn(t)l\\X„{t)\\ when ||X„(f)|| > 0, 
and also Ln{t) - n~^ ( ni.n^ t) . In the sequel, we denote by iJ a Brownian motion with 
drift -Aa, variance 2A; P'' for £ IR denotes its law started at b. Note that it is known 
that under the heavy- traffic assumption, the sequence (L„) converges weakly to B_. The 
goal of this section is to prove the following result. 
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Theorem 5.1. The sequence of processes {Xn, n>l) underP^ converges weakly as n goes 
to infinity toBn under P^. 

The above theorem is reminiscent of heavy- traffic diffusion limits for the joint queue 
length process in various queueing networks, involving reflected Brownian motion and 
state space collapse, see for instance Bramson 1 3 1 , Reiman 1 13J , Stolyar [16} and Williams il9J . 
However, to the best of our knowledge, this is the first result which shows that mobility 
of users, rather than scheduling, routing or load balancing, can act as a mechanism pro- 
ducing state space collapse. 

It is straightforward to adapt the proof ofTheorem lS.ll to handle a more general initial 
condition. Let b>0 and assume that (0) || — ► b, then it can be proved that: 

• for any sequence such that e„ > 0, e„ ^ and ?z^e„ +oo, the sequence of 
shifted processes Xn) converges weakly towards Bn under P^; 

• if Rn (0) 71 then (X„] converges weakly to Bn under P^. 

We see that if i?„ (0) n' ^ n then converges in the sense of finite-dimensional 
distributions to a process which is discontinuous at and so cannot converge weakly 
(at least in the space of cadlag functions). 

5.1. Overviewof the proof. Remember that e\{Xn) - i<^°0THx„.e)^^'^ri) and that geiX,^) 
is the left endpoint of the first excursion of Xn with height larger than £, equivalently 
the left endpoint of the excursion of X„ straddling r' To prove Theorem l5.ll we 
use Theorem 4 in Lambert and Simatos |8|: this theorem states that the conclusion of 
Theorem lS.ll will hold if we can show that g^iXn) => ge{Bn), e\[Xn) => e\[Bn) and (To o 
s\) [Xn) => ( To o el) [Bn) for any £ > 0, where X„ is considered under P° and B under P"; 
note thai gir[Bn) = g£.(_B) and similarly that [TQoel)[Bn) = [TQoel)[E). 

The convergence of the two sequences [e\[Xn)) and [[Toe\)[Xn)) is studied inPropo- 
sition l574l its proof relies on the following ideas. First, let G)n - Xn[T^ [Xn,cn)): the 
Markov property shows that e\[Xn) under P° is equal in distribution to (7(X„) under 
P^". In Lemma [5^ we show that a),j =; eu, i.e., the system is with high probability ho- 
mogenized at time r' [Xn,En). Then, we want to show that this property lasts on times 
of order of n^ that we are interested in. To do so we exploit the bound of Proposition l4.2l 
which shows that we wUl be able to control the process on times of order of n^ as long as 
there are at least of the order of log n users in the network. Thus, this reduces the prob- 
lem to control a[Xn) started with logn users: this is studied in Lemma [53l where it is 
shown that ro(X„) is at most of order of Since arrivals and departures are Poisson, 
on this time scale the total number of users at most varies by \/n which is neglegible 
compared to the space scale n that we are interested in. 

On the other hand, the convergence ge[Xn) ^ ge[Bn), proved in Proposition l5.5l fol- 
lows from two arguments: one that provides a lower bound in terms of a sum of i.i.d. 
terms related to (To o e\)[Xn), whose asymptotic behavior we will control thanks to the 
convergences of [el[Xn)) and [[Tqo el)[Xn)); and one that provides a corresponding up- 
per bound via the coupling with the Ml Ml 1 queue. 

5.2. Convergence of (gglXn), « > I) and ((Too e|.)(Jf„), « > I). The convergence of these 
two sequences is proved in Proposition l5.4l 

Lemma 5.2. Forany£,S>0, 

lim Pl[\\R„[ThXn,£))-n\\>6]^Q. 

n— >+oo V ; 
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Proof. Fix c > 0, by monotonicity, the result only needs to be proved for 5 small enough; 
in the sequel we will consider 5 > such that rj < 2eo with rj - nS^ISZ. Let t - TirjUAK)), 
Un - n^'^, Vn- n^'^, Tn - [Xn, tie), T'^ = r' [Xn, ne- Un) and £„ be the event 

£„ = { r; + T < r„} n { II x„ (r; + T) II > - 2m„} n { II r„(r; + T) - ;r|| < 77} . 

Using that T,', < Tn, the strong Markov property applied at T'^ gives 

V>l[E'„)^Vn"[{Tn<T}li{\\Xn{T)\\<nE-2Un}li{\\rn{T)-n\\^r]}) 

where (D'„ is equal in distribution to Xn{T'^) under P°. On the other hand, the strong 
Markov property applied at the stopping time T'^ + t gives 

P° (III VniTn) - n\\ > 5} n En) < max P^, (|| r„(r„) - ;r|| > 5) . 

with 3'n - {y e : ne - 2m„ < ||y|| < ne and g{y) < 77}. Let y„ e STn that realizes the 
maximum, so that P° ({|| r„{T„) - 7r|| > 5} n £„) < P^" (|| r„(r„) - 7r|| > 5). Thus we get the 
bound 

P° (II ( r„) - ;r II > 5) < P^" ( r„ < T or II x„ (T) II < - 2 M„ ) + pf" ( || r„ (t) - ;r || > 77) 

+ P^"(ll'-„(r„)-;r||><5). 

Since ||(D'„|| = \n£- Un], if under P^" we have r„ < t or ||x„(t)|| < ne-2un then we 
must have at least m„ arrivals or m„ departures in [0,t]. Since the arrival and depar- 
ture rates are bounded (by jc) while m„ +00 we see that the probability of this event 

vanishes, i.e., P5?"(r„ < t or ||x„(t)|| < ne-lun) 0. 
On the other hand, since ||(D'„|| = \ne- m„1, wehave 

P^" (II Tn (T) - 71 II > 77) < P^" (II Tn (t) - 71 || > 77) 

for some yj, e with ||yjj|| - \ne- Un]. Since 77 < 2eo and t - TirjliAK)), Lemma l4Jl 
shows that Pj^" (|| r,, (t) - :7r|| > 77) ^ 0. Hence we have 

limsup po (II r„(r„) - ;r|| > 5) < limsup P^« (|| r„(r„) - ;r|| > 5) 

and we now show that this last upper bound is equal to 0. 
We have 

(II rniTn) - n\\ > 5) < P^" (r„ > y„) + P^" (r' (r„ - 7r,5) < 

and by Lemma |2!T] 

P^" (T„ > Vn) < P^" (r' {(n, riE) > y„) < Pf (r' (^„, nc) > 

where yJJ e is such that ||yj,'|| - [ne- u„J, using stochastic monotonicity of {(n, ne) 
in the size of the initial condition. Since £ fi ^ i we get 

P^/' {Tn > y„) < V>i" [t^ (In, ne) > = pO (r' u„) > y„) . 
On the other hand we have 
Pi" [Thrn-n,S)< y„) < P^« ( T^ [In, Un) < + P^" ( r' (r„ - ;r, 5) < A r' (x„, 

where we have used the inequality (x„, < t' m„) that stems from Lemma [2!T] 
Since 77 < Icq, Q[yn) £ f] and ||y„|| > Un, all the assumptions of Proposition FOl are sat- 
isfied, at least for n large enough. Thus the second term of the right-hand side of the 
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previous display vanishes, and gathering all the previous bounds we see that we are left 
with 

limsup P° (||r„(r„) - n\\ >8)< limsup P° (r' u„) > Vn] 

n— '+00 n^+oo ^ ^ 

+ limsup P^" ( r' u„) < J. 

Since the sequence of rescaled processes (L„) with L,j(f) - (!n{n^t)ln converges in 
distribution to a Brownian motion, it is not hard to show that the two sequences of ran- 
dom variables [Un^T^ [in, Un)) Under P° and [rT^T^ {£„, u„]) under P^" converge weakly 
to a non- degenerate random variable (actually, hitting times are continuous function- 
als when the limiting process is the almost sure realization of a Brownian motion, see 
for instance Proposition VI. 2. 11 in Jacod and Shiryaev [7]). Since « Vn « this 
finally proves that the right-hand side of the previous display is equal to 0, hence the 
result. □ 

Lemma 5.3. Letcpn - [(\ogn)^\: then 

lim max P>UTo[x„) > y/n] \ ^0. 

n-+oo^y:||y||=<^„ 'j 

Proof. Let M„ £ N be such that 2^""^ <(t>n^ 2^" ■ By monotonicity of Tq in the size of 
the initial state, we have 

max P^(ro>v^< max Pi(To>\^ 

y:\\y\\='Pn y^3-[M„) ' 

where from now on 3'{m) = {y e : ||y|| - 2'"} and we omit the dependency of the 
functional operators when they are applied at x„, so that To - ro(x„). Define S{f) - 
if, 11/(0) II /2) and (pn,m = exp(M„ + 2'-""-'"^"^], and note that since p„ < 1 by assump- 
tion, S is almost surely finite since x„ is recurrent by Proposition 14.31 The relation 
Tq - S + TooQg gives for any n>l and y e 

Pi (ro > v/?i) = Pi {To >VJi,S^ (pn,o) + Pi {To >VJi,S< (pn,o) 

<Pi{S> (pn,o) + P^ (To o 0s > V^-cpn,o) ■ 

Since moreover Pj^(x„(S) £.9" (M„-l)) - 1 for y £ 5" (M„), the strong Markov property 
at S gives 

max Pi{To>Vn]< max Pi{S>(p„,o)+ max Pi{To> Vn- (pn.o) ■ 
Iterating M„ -i- 1 times, we obtain 

max P^, (To > < X max P^^ (S > q)n,m) 

+ niax Pl{To>s/n-(pn,o (Pn,M„]- 

y^<J (0) 

Now 5^(0) = {0} and To under P° is equal to 0. Since in addition 2*^" < 2(pn £ 2(log?i)2, 
we obtain 1^""^ <2^"^ (log n) ^'^ and so 

<Pn,o + • • • + (Pn.M,, ^ i-Mn + 1) exp (M„ + 2^"'*) < exp ((1/2) log n)^^/n 

for n large enough. For those n, we get 

max Pi{To>^/n-(pn,o <Pn,M„)=0 

ye3"(0) 
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and finally, we obtain after a change of variables 

M„ 

(9) max IP^(ro>^<y max Pf, (S > exp(M„ + 2'"'^)] . 

ye3-(M„) yeS^C") 

Let us justify that for fixed m > and y e we have 

(10) n^^^oo"^" ~ '^^P^^" + 2"""^)) = 0. 

First, let {u{n)) be a subsequence such that 

lim [S > exp(M„(„) + 2""^)) = limsup (S > exp(M„ + 2'"'*)) . 

Since the [0,cxd)^^ -valued sequence {(.^uin),k< l^u(n),ky 1 £ £ « > 1) lives in a com- 
pact set as a consequence of the heavy-traffic assumption, we can find a subsequence 
(.v{n)) of (w(?z)) and Aoo-Moo e [0,oo)^ such that fc ^ Aoo,fc and Ai„(„),fc ^ f^oo.fc for 
each A; = 1,...,K. Because of the heavy- traffic assumption, we have || Aqq II = II Moo II ■ 

It is then not hard to see that the sequence (Xj/f/^)) under converges weakly to Xqq, 
where Xoo is the Markov process with Xoo (0) - y and generator D.ao defined similarly as 
Q„ but with A„jt and idn,k replaced by Aoo.k and fioo„fc. respectively. Since Xoo lives in 
and is piecewise constant, it is not hard to prove that S is a continuous functional 
at Xoo and so the continuous-mapping theorem implies the weak convergence of the 
sequence (S(Xy(n))) towards S(Xoo)- In particular, 

- '^^P^*^"'") + 2'"'^)) = P (S(Xoo) = +00) . 



Since ||Aooll = IIa*ooII. Xqo is recurrent by Proposition fOl and so S(Xoo) is finite almost 
surely which proves fTO) . Since for each fixed m > 1 the set ^ ( m) is finite, combining (9) 
and (To) we obtain for any M > 

I \ ( M„ \ 

limsup max P]^ ( To > < limsup Y max P^ (S > exp(2'"'*)) < Y Um 

n^+oo Vye3-(M„) j n^+oo \m=M>'^-^''"' / m>M 

where 



{/,„ = sup[ max P^(S>exp(2'"'^))l. 

n>l Xy^-Tim) j 



Thus if we can prove that the series [Um) is summable, letting M +oo in the previ- 
ous inequality will show the result. 

For the rest of the proof, fix any 5 > such that S <n and rj < 2eo where rj - nS^ 132. 
Let in addition t - TirjliAK)). Then for any y e N^, we have 

Pi [S > exp(2'"'*)) <Pi{S> exp(2"^'*), || r„ (t) - ;r || < 77, | || x„ (t) || - || x„ (0) || | < m) 

+ P« (II r„ (T) - ;r II > 77) + P^ (I II x„ (T) II - || x„ (0) || | > m) . 

Invoking the Markov property at time t for the first term and l|4) together with Lemma lZ2 
for the two last ones, we get 

max P^ (S > exp(2'"'*)] < max P^ (S > exp(2'"''') - t] 

ye^im) ^ ' ye3-'{m,T]) ^ ' 

+ max P^f||r„(T)-:^|| >77]-i-P(3»(k:t)> m) 

y£.T{m) 

where 5''{m,8) - {y eN^ : \\\y\\ -2"'\ < m and g{y] < 77}. The last term of the above 
upper bound defines a summable series since a Poisson random variable has a finite 
mean; the second term also by Lemma fOI since 77 < 2£o. It remains to control the first 
term. 
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So consider y e Sr'{m,S) and let Vm - expCZ*"'*) - t and ^Jm -2^ ^ - mil. Since 
llyll >2'"- m wehave S= rUx„, |ly|i/2) < THx„,iffm) under Pj^ and we obtain 

P^ (S > Vm) < P^ ( U„, i/^m) >Vm,Thr„-7i,S]>vJj 

+ py[Thxn,Vm]>Vm,Thrn-7T,5]<Vni). 

Since 5 < tt, we have ||x„(r)|| - (nit) for all f < T' (r„ - :fr,5) by (3) and so 

V>i[Thxn,Vm)^Vrn,Thrn- 71,8)^ Vrn]^K[Th{n,Vm)^Vrn,Thrn- 71,8)^ Vrn] 

which implies 

max ¥'i{S>v,n)< max P^ (rU/„,i/A^) > y J 

y£3~'(m,i]) ye5"'(m,»)) ^ ^ 

+ max py(Thrn-7T,S)<VmAThxn,V/m)]- 

y£3~'{m,r]) 

As for the second term, it is easily checked that all the assumptions of Proposition l4.2l 
are satisfied, at least for m large enough, since 8 and 77 < 2eo are fixed and both Vm 
and 2"^"^ grow without bounds with m. Proposition 14.21 provides a bound uniform in 
y e 3~'[m,8) and n> \ which defines a summable series in m. Hence to complete the 
proof it remains to show that 

(11) ^ SUp{p^"'(rk^„,Vm)>!^m)}<+CX) 

m>Mn>l 

with ym - (2'" + m)ei, using monotonicity of i(n,Vm) in the size of the initial state 
and the fact that £n only depends on the initial state through its size. Note that 

P^ [Th(n,Vm) > Vm] = P^ [ThIn,Vm) > ^m) = P° (rk/„,Vm) ^ 

with = - 2"" - m = 2'""! - 2"" - 3m/2. It is well-known that TH?n,f'm) un- 
der P° scales like (v'm)^ ~ 2^"* for large m and n, which is negligible compared to Vm ~ 
exp(2'"'4). In the case pn - 1, we can use the reflection principle and establish an expo- 
nential upper bound as we have done in the proof of Proposition l4.3l When Pn < 1 we 
can couple ^„ with a critical random walk £'„ such that > In particular we have 

THfn,Vr„)<THf'n,Vm) and SO P^'" (fl (/„,Vm) > !^m) < (^Vm) > !^m), where 

this last term obeys to an exponential upper bound. This proves fTTi and compeletes 
the proof. □ 

Proposition 5.4. For any e > 0, the sequences [eliXn), n> 1) and [{Tq o e\){Xn), n> \) 
underP^ converge weakly to el [Bn) andiTQoel){B) underP^, respectively. 

Proof. In the rest of the proof fix any e > and let a)„ be the law of x„(r' (x„,£)) un- 
der P° . In view of the strong Markov property, the convergence properties claimed are 
equivalent to the convergence of the two sequences {a{X„),n > 1) and (ro(X„), n> I) 
under P^" towards oiBn) and To(B) under P*^, respectively. 

We know that ((T(L„:/r)) and {ToiLn)) under P^" converge weakly to a{Bn) and Tq{B) 
under P*^, respectively, and we want to transfer this result to X„ using Theorem 3.1 in 
Billingsley [I], sometimes referred to as a "convergence-together" result. Thus we only 
have to prove that for any ;0 > 0, 
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Define (/)„ - L(log?z)^J as well as Tn,x - iXn, n ^(pn) and Tn,L - r' (L„, n ^(pn)'- then 

mXn) = Tn,X + [To o9T„_^)[Xn) and ToiLn) = Tn,L + [To oBT„,)[Ln) 

so that 

K" (I To[Xn) - To[Ln)\ >P)< Pn" {\T„,x " T„,l\ > O) 

+ Pt" {[To ° 0r„,x)(^«) 5 i3/3) + P^" {[To °OT„^,)(Ln) > iS/3) . 

Note that Lemma l2Jlimplies that P^"(|r„,x- Tn,L\ > 0) = P^"(r„,x > r„,i). Together 
with the strong Markov property at Tn,x and Tn,L, this gives 

K"{\TQ[Xn)-To(Ln)\>P]<Pn"{Tn.x>Tn.L]+ max Pf, (ro(x„) > n^/j/a) 

y:llyll=</>n 

+ max P^(ro(^„)>«2/3/3). 

y:||yll=</>n 

Since the sequence [(p~^ To(£„)] under P^J""' converges in distribution to a non-degenerate 
random variable, the last term goes to since » cj)^. The second term goes to by 
Lemma [531 and so it remains to control the first term. In the rest of the proof, let 5 > 
such that d<7£ and rj - nS^ 132 < 2co : we have 

Pn" {Tn.x > Tn,L] < ^t" (ll^n(O) - ;r || > r?) + max P]; {Tn,x > Tn,L] 

y^STn 

with STn = {y e : ||y|| = \_ne\ and Q[y) < 77}. The first term goes to by Lemma W2\ and 
so we have at that point 

limsup P^" (|ro(X„) - ^o[Ln^\ >P]< llmsup (max P^(r„,x > Tn.L]] ■ 
+00 +00 Vye3"„ ; 

Thanks to Lemma lZTl one sees that Tn,x > Tn,L implies that there was a time t < Tn,L 
such that Xn,k[t] - for some A;. At that time we have \\rn{t)-n\\ >n>S and so 

Pi {Tn.X > Tn,L) £ (r' (r„ -7l,5)< r,,,^) . 

Further we have 

pI{t^ [rn -7i,d)< T„,l] <Pl{n< Tn.i] + Pn {t^ [r,, -n,6)<nA T^.l) . 

The first term goes to uniformly in y with ||y|| = YEn\ since the sequence [Tnx) 
under P^ with y^STn converges in distribution, while the second term can be rewritten 
as 

Pi[Thrn-n,n)<nf\ Tn,i] = P^ ( r' (r„ - ;r, 5) < a rk^„, 

which goes to uniformly in y e 3'n by Proposition |4j2] (using r' [£n:<Pn} ^ T^ [Xn,(pn))- 
This proves the result on ToiXn), the result on cr(X„) follows along the same lines but 
at the expense of more technical details. Lemma [ZTI and (3) imply that Tn,x - Tn,L and 
Ln[T„,L) = \\X„[Tn.xn ^n-^(p fi in the event {T^ [Rn — 71,6) > Tn^i], so that the strong 
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Markov property gives 

Pt" [sup\\a{X„){t)-(j{Ln7im\\ > js] < Pt" [ThRn-^,S) > 



sup WXnit)- Ln{t)n\\> p/3\+ max Pi 



0< t< T„ 



y£<S>„ 



sup ||X„(f)||>/3/3 

\o<t<ro(Jf„) / 



+ max I 



sup L„(f)>)0/3 
\o<i<roa„) / 



where 0„ = {y e : ||y|| - (pn). We have already proved earlier in the proof that the first 
term vanishes. The last term can be seen to vanish invoking the convergence of (L„) 
towards a Brownian motion. The third term vanishes because for any y e <!>„, 



sup ||X„(r)||>;6/3 

\o<t<ro(x„) 



sup ||x„(r)|| > 
\o<t<roU„) 

< Pi [MXn) > + > «iS/3 - cPn] 



and both terms go to uniformly in y £ $„, the first one using Lemma [531 and the sec- 
ond one using Markov inequality. It remains to show that the second term also vanishes. 
Let Sn - sup[Q j.^^] L„: then for any < t< Tn,L, one has 

II X„ ( - L„ ( r) ;r II < II X„ ( II II i?„ ( r) - TT II + II ;r II I II X„ ( II - L„ ( r) I 

< L„(f)||J?„(r)-;r||+2|||X„(r)||-L„(f)| 

< S„||i?„(r)-;r||+2|||Jf„(r)||-L„(f)| 

where (1) follows by adding and subtracting n^||X„(f)|| and using the triangular inequal- 
ity, (2) follows from || X„ ( f ) || < L„ ( f) -i- 1 1| X„ ( r) || - L„ { f ) | together with || i?„ ( f ) - :/r || < 1 and 
\\n\\ - 1 and (3) is by definition of Sn, since t < Tn,L- Using this upper bound together 
with standard manipulations, we get the following upper bound, valid for any s > 0: 



sup ||X„(f)-;rL„(f)|l>/3/3 

o<t<r„,i 



sup \\Rn[t)-n\\>p!{&s] 



sup |||X„(r)||-L„(f)|>0 

o<t<r„,L 



Similarly as before the first and last terms go to 0, so we are left with 



limsupP^" 



sup \\Xn{t)-nLn{t)\\^p 



<limsupP®"(S„>s). 

n— >+oo 



Letting s +oo completes the proof, since the sequence (S„) under P^" converges 



weakly to sup 



|o,ro(B)i 



B under P*" 



□ 



5.3. Convergence of {gi;iXn),n > 1). To complete the proof of Theorem 15. 1 1 based on 
Theorem 4 in Lambert and Simatos |8| it remains to be shown that gdXn) ^ geiM)- 

Proposition 5.5. Forany£> 0, the sequence [geiXn)) underP'^ converges weakly to gs(B) 
under P°. 
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Proof. Since > L„ by Lemma [ZTl it is clear that rUx„,£) < T^[Ln,£), and hence 
going back in time and using again || X„ || > L„ we see that ge [Xn] < geiLn). On the other 
hand, since L„ => B it is not difficult to show that ge{L„) => gelB), see for instance Lam- 
bert and Simatos [9J where similar computations are carried out. This proves that the 
sequence igeiXn)] is tight and that any accumulation point is stochastically dominated 
by ge{B). We now derive a corresponding lower bound which will conclude the proof. 

Let d > and denote by fij^ ^ the fcth excursion of Xn that reaches level S shifted at 
the first time it reaches this value: formally, we have 



X- sUXn) and e\ , „ : 

n,l,5 "■' ;j,fc+l,o 



[^l°ea„,,,]{Xn) 



with dn,fc,5 the right endpoint of the excursion of Xn corresponding to En,k,s- Let Nfj^s.E 
be the number of excursions of X„ that reach level S and not level £ before the first 
excursion of X„ to reach level e: it satisfies 

A^„,5,£ + 1 = infi fc> 1 : supllfi), „(f)|i > c 

Then Nfi,s,e is a geometric random variable with parameter Pn,s,£ given by Pn,s,£ - 
P°(sup||eg(X„)|| < £), the ^.^,1 < fc< N„^s,e) are i.i.d., independent ofNn,s,£ and with 
common distribution eg conditioned on {sup || || < e}, and we have 

gAXn)> E To <J 
with the convention - 0. This last inequality implies for any s > 

and it can be computed that this last upper bound is equal to 

l-Pn.S,e P°(sup||4(^„)ll>£) 



1 - pn,s.eEl (e-^<^°°''«»^"' I sup||4(XJ|| <e] 1 - E" (e-^'^°°''«»^"';sup||4(X„)|| < e] ' 

Proposition l5.4l and the continuous-mapping theorem give 

^liin^P'" (sup II 4 (Xn) II >e]^P° (sup4(:B) > e] . 

On the other hand, the two convergences of Proposition 15.41 can be shown to hold 
jointly, see for instance Lambert and Simatos i8J, and so the continuous-mapping the- 
orem gives 

Thus for any 6> we have 

P° (sup II 4(B) II 5 £) 



^lirn^E° (e-^'^°°4)(^");sup||4(X„)|| <e]^E° (e-^'^o^P^jsup ^(B) < e] . 
Are have 

limsupE°(e-'«^t^''') = 



.Eo(e-^(ro°4'®;sup||4(B)||<£) 

irsion theory, the abov 
md so we finally get 

limsup E° (e-^S' '^"') < E° (e"^^-® ) . 



Using standard arguments from excursion theory, the above upper bound is seen to 
converge towards E''(e~''^^'-') as 5 ^ 0, and so we finally get 
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This implies that any accumulation point of the tight sequence [gsiXn)) is stochasti- 
cally lower bounded by geiB), and since a corresponding stochastic upper bound holds 
this gives the result. □ 



6. Convergence of the stationary distributions 

Throughout this section the heavy- traffic assumption of Section[5]continues to be in 
force, and we moreover assume that a > 0, so that Pn < 1 for n large enough. Recall 
that in this case x,, is positive-recurrent, see [51115], and that v,j denotes its stationary 
distribution. This section is devoted to proving the following result, where we write 
similarly as in the previous section Xf,{0) - «~^x„(0). 

Theorem 6.1. The sequence (X„(0), n> I) under P)^" converges weakly as n goes to in- 
finity to En where E is an exponential random variable with parameter a, and all higher 
moments converge as well, i.e., Ej^" (0) || — ► r\la^ for all integer r > 0. 

Since the exponential random variable E in the above theorem has the stationary 
distribution of the reflected Brownian motion B_ introduced in Theorem l5.1l we may ob- 
serve that the heavy- traffic characteristics are preserved under an interchange of limits. 
While such an interchange of limits tends to apply in most specific cases, there do not 
appear to be any general guarantees for that. The proof of Theorem 16.11 relies on the 
following estimate. 

Lemma 6.2. There exist two constants c, d £ (0,oo) such that for every n>l, every q>0 
and every l< k< K, 

Pi" (II x„ (0)11 > q,Xn,km = 0) < ce-''l. 

Proof. Fix n > 1, q > and 1 < k < K. Recall the constant eo of Lemma ISTI and let 
< £ < 1 be any number such that 

1 < < 1. 

n 7r(l-e) 

Let 77 = 2£/(7r(l-e)), so that < ?? < 1 and (,l-r])n<Eo, and t = T{[l-r])n/{2K)). Then 
for any y e N^, one has 

(lU„(T)|| > q,Xn,kW^O)<Pi{\\XniT)\\ > ^7, r„,fc(T) < £71^) 

< (II Xn (t) II > q, Xn,k (t) < C^Tfc || x„ (t) || , an (t) V dn (t) <Eq) 

+ Pi[an{T)\J dn{T)>£q]. 

The last term P^(a„(T) v (i„(T) > eq] is upper bounded by P(5*(7ct) > eq) which we 
control using (1). One can check that h{x) > x for x > e, hence for q such that exi < eq 
we obtain 

Pn(«nW vd„(T) >Eq] < e~^'^ , q> CKT/e. 

In particular there exists a finite constant c such that P^(fln(T) v dnM > eq) < ce"'^'? 
for all q > 0. On the other hand, {||x„(t)|| > q.anM < eq] c {||x„(0)|| > (1 - e)q} and 
according to ID, we also have the inclusion 

{fl„(T) vd„(T) < eq} c |x„,fc(T) > x^_j.(t) -£^,||x„(t)|| < ||x„(0)|| -H£^|, 
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hence 

( II Xn (t) II > q, x„^k W ^£^k\\ Xn (t) II , a„ (t) V d„ (t) <£q] 

<Pi[\\ Xn (0) II > (1 - £) q, J. (T) <E7lk{\\ X„ (0) || + + 

< (ll Xn (0) II > (1 - E) q, r;_j. (T) < e;rfc (1 + £/ (1 - £)) + £■/ (1 - £)) . 

Recalling that 77 = 2e/ (tt (1 - e)) , one sees that the last term of the previous equation is 
upper bounded by l{|iy||>(i-E)(/}P^(r^ j.(t) < rjjik) and we have further 

l{\\y\l^a-e)q]Pi {r'n,kM < V^k) £ l{iiyii>(i-£)<?![P^ (ll r'nM - ;r|| > (1 - 77)71) . 
Since < (1 - ri)n < eq, Lemma lSTTI implies 

h\\y\\m-e)q}Pl (II r'nM " ;r|| > (1 - Tfin) < 2^:exp [ '-jj^^ . 

This proves the result, with for instance c' = (2K]~^mm[4K^E, (1 - rji^jf'il - £]) and c - 
C + 2K. □ 



Proof of Theorem \6J[ We first prove the convergence of the moments by induction on 
r > 0. The result is immediate for r = so consider r > 1 and assume by induction 
hypothesis that E^" (||X,j(0)|p) — ► s\/a^ for every < s < r - 1. Let G„ be a geometrically 
distributed random variable with parameter p„: Lemma lOl implies that 

E;;"(iix„(o)r)>E(G;,), 

and using 77(l-p„) ^ a, it can be proved that EUGnlny) —•■ rl/a'^. This provides a lower 
bound and so we only have to show that limsup„^+j^E^"{||X„{0)ir) < rlla''. Let m> 1: 
summing the balance equations over the set {y e fM^ : || y || < m - 1} yields 

(12) A„C (II x„ (0)11 = m - 1) = i^nK" illxnion = m) 

K 

- Z ^'n.kPn' (lUn(O)ll = m,X„,fc(0) = O) . 
fc=l 

For any M>1, writing m'' - Z^^o siJ-syS '^ ~ 8^^^^ 

M r y\ 

X '"''C (lUn(O)ll = m- 1) = ^ ' ,, E;" (||x„(0)r;||x„(0)|| < M- 1) 

and so multiplying {12) with n/ 1 iin on both sides and summing over 1 < m < M, we 
obtain 

Z ,r .. PnE^" (||x„(0)r; IU„(0)|| < M- 1) = e;- (||x„(0)r; |U„(0)|| <M) 

s=o Sif- •?)! 

1 ^ 

E /"n.fcE;;" (||x„(0)ir;||x„(0)|| < M,x„,fc(0) = 0). 

fc=l 

Using e;;" (||x„(0)r;llx„(0)|| < M) > e;;" (I1x„(0)||''; ||x„(0)|| < M- 1), other simple in- 
equalities and isolating the terms corresponding to 5 = r and 5 = r - 1 in the previous 
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sum, we end up with 

r-2 

{l-pnWl" (ll^n(O)r; llxn(0)|| < M- l) < rE^" (||x„ (0) || + r! ^ E^" (||x„(0)f ) 

s=0 

K 

+ Y.El"[\\xnmV;xn,km^o]. 

k=l 

Letting M +00 and dividing by n''{l - p,,) gives 

(ll^«(0)r) < — -El" {\\Xnm\r') + —7^—— E "'"""""e;" (iix„(o)f ) 

n{l-pn) na-pn)f-^o 

1 ^ 

+ ^7; 7 E En" (II ^n (0) ir; (0) = O) . 

«'(1-Pn) fc=i 

Since nil - p„) — ► a > 0, by induction hypothesis the first term of the above upper 
bound converges to r!/a;' while the second term vanishes, thanks to the terms n*"'''"^' 
that go to for 5 < r - 2. We now show that the last term also vanishes, thus concluding 
the proof. To this end we prove that Ej^" (||jc„(0)ir;x„,fc(0) = O) ^ for all 1 < k<K. Let 
ke {!,..., K] and M > 0: we distinsuish the two events {||jc„(0)|| < M] and {||x„(0)|| > 
M + 1}. On the one hand we have 

E^" [Wxnmf-AxnmW < M,X„,fc(0) = 0) < M'Pl" (II x„ (0)11 < M) < M'' [l-[pn)'^) 

invoking Lemma lOl to get the last bound. On the other hand, Lemma lE2l gives 

(lU^COir; ||x„(0]|| > M+ l,x„,fc(0) = 0) = E ^n' (lUn(O)r > m,x„,fc(0) = O) 

m>M+l 

< E ce-''""'"'. 

Since pn^l, letting « ^ +00 for a given fixed M gives 

limsupE;,"(||x„(0)||'';x„,fc(0) = 0)< E ce"'^''""'- 

Letting M ^ +00 achieves to prove that integer moments converge. 

We now prove the weak convergence result. Since Ej^"(||Jf,3(0)|p) 21a?-, the se- 
quence (||X„(0)||, « > 1) is uniformly integrable and tight. Let X be any accumulation 
point and assume without loss of generality that Jf„ (0) ^ X: we wiU prove that X - En, 
thus proving the desired result. 

Since n{\ - pn) a, \\X\\ is stochastically lower bounded by an exponential random 
variable with parameter a thanks to Lemma [2J] On the other hand, since X„(0) ^ X, 
En" (II -^n (0)11) lla and (||X„(0)||, n > 1) is uniformly integrable, we get convergence of 
the means and so E(||X||) = l/a. In summary, if E is an exponential random variable 
with parameter a, then ||X|| is stochastically lower bounded by E and E(||X||) - E{E): 
hence ||X|| and E must be equal in distribution. 

In particular, P(||X|| = 0) = and so the continuous-mapping theorem implies that 
||i?„(0)-7r|| ^ ||X/||X|| -:^||. LetO<e<2£:o and t = T(e/(4i<:)): since 

K" iWRnm - ^11 s £) = Pi" (II r„(T) - ;r|| > £) = 

Lemmal4Jl implies that Pi" (||i?„(0) -7i\\>e)^0 and so P(||X/||X|| - 7r|| > e) = 0. Since 
£ < 2£q is arbitrary, letting £ ^ shows that X = ||X||:^ and since ||X|| is an exponential 
random variable vnth parameter a this proves the result. □ 
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7. Convergence of the sojourn times 

Throughout this section the heavy- traffic assumption of Section |5] continues to be 
in force, and we moreover assume that a > so that Pn < 1 for n large enough. We 
fix a sequence (y,,) in with p(y„) and ||y„||/n be (0,oo), and we return to the 
network description of x,, , see Section [Z2l In the sequel we implicitly consider x„ under 
Pn" and write => for weak convergence. 

For n> 1, we pick an initial user of the nth system uniformly at random, which we 
refer to as the tagged user, and denote by £„ its initial service requirement (which is 
exponentially distributed with unit mean), by^^„ its trajectory and by x„ its sojourn time, 
so the following relations hold: 

X„^mf{t>0:sn[t)^ En} with s„(f)= f - — du. 

Jo 1 VX„,f„(„)(M) 

Note that 5„(f) for r < Xn is the service received by the tagged user up to time t. We 
want to show that rT^Xn ^ bEI X where in this section E is an exponential random 
variable with unit mean. Introducing the rescaled service process S„(t) - Snint), we 
have n~^Xn = inf{f > : Snit) - £„}. Let s'„ be defined as follows: 



s„(f) = - — - / di 

II yn II JO 7rf„(„) 



and define S'„ {t)^s'n{nt). Let finaUy S = (fo" U f , f > 0) . 

Lemma7.1. The sequence of processes (S',^) underP^" converges weakly to S . 

Proof. Since S'„ and S are strictly increasing and continuous, we only have to show that 
finite-dimensional distributions converge, see for instance Jacod and Shiryaev |T The- 
orem VI. 3. 37]. Since the law of ^„ does not depend on n, we can couple all the processes 
(S'„, « > 1) on a common probability space by taking ^n- - Then we have 

llynlUo T^iiiu) II yn II it=l T^k JO 



and so 



s'„(f)- 



~b 



1 



llynll n 



71 1 Jo 



,lc}du-n^nt 



llynll 



The second term goes to by assumption while the first term is easily seen to vanish 
using the ergodic theorem. This shows that S'„{t) S(.t) almost surely, for every f > 0, 
which readily implies convergence of the finite-dimensional distributions. This proves 
the result. □ 



Lemma 7.2. The sequence of processes (S,,) under converges weakly to S. 

Proof. Since S',, ^ S by Lemma [7J] Theorem 3.1 in Billingsley [IJ shows that it suffices 
to prove that for any e > and any tQ>0, 

Urn ¥>l"\ sup |s„(f)-s'(f)|>£ =0. 

n-+oo \.o<t<nto ; 
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By definition we have 

\Sn{t]-s\ 



\J0 JCn,^„(i()Vl llynllJo ^(„iu) 



< 



1 P-n,i„{u) I II III J 

ynW Jo yif„(u)iXnI„iu)^l) 



WynW Jo 7r^„MU«,f„(u) Vl) 

Fix in the sequel c = 2 + sup „> j ( II yn II / «) and 5>0 such that 

„ „ I l^nScn^to \ 
S<n, 8<bnlc, < 2£o and sup \<£l2, 

~ ~ n7il\\\yn\\E^E\\yn\\-0cn)j 

where rj - n8^l?,2. Let F„ be the event 

= { r' (r„ - TT, <5) > « fo } n { r' (II x„ II - II y„ II , n5) > n fo } 

where ||jc„|| - ||y„|| refers to the process (||x„(f)|| - Hynll, f 5 0). Since S < n we have 
r'Cn - ^,S] < To{x„] by (3) and so x„_fc(f) > for all t < ntQ and 1 < /c < in P„, in 
particular Xn,k{t] v 1 = Xn_fc(f). Thus, we also have in F„ 

|x„_fc(f) vl-:?rfc||y„||| = |x„,fc(f) - n-fc||y„|| | < |x„_s;(f) -7rs;||x„(f)||| + 7rfc |||x„(f)|| - ||y„||| 

s (llynll + 5n) ||rn(t) - ttII + 5« < 5c« 

for all t < nto and 1 < fc < ^C. In particular, x„ j;(f) > 7r||y„|| - Sen which is positive for n 
large enough since dc< bn. Plugging in these different bounds, we obtain the following 
inequality, valid in F„ for any f < nto: 

I ; I 1 r'^ Hn „ Linden^ to 

\sn{t)-s'{t)\< / — 8cndu< — <£/2. 

Ilynll Jo 7r(7r||y„|| -dc«) llynllzr(Z£lly«ll - ocn) 

Thus we obtain | s„ ( f ) - sj, ( f ) | < c in f „ and therefore 

Vl" [ sup 1 5„ ( r) - ( f) I > e] = (F^) < P^" ( r' (II x„ II - II yn II ,n8)<nto] 

Vo<t<nto ' 

+ P^"(r'(r„-;r,5)< nfoAr'(llx„||-||y„||,«5)). 

Since T' ( II x„ || - || y„ || , n8] < (x„ , || y„ || - n8) , the second term can be shown to go to 
by Proposition 14.21 fusing that giyn) < rj ioi n large enough since giyn) 0). The first 
can also be shown to go to 0, since 

P^" (r' (||x„|| - ||y„||, ?25) <nto]^ P^" (r'(IIX„|| - ||X„(0)||,<5) < to/n ) 

and the convergence Jf„ ^ Bji of Theorem [STTI implies that r'dl^nll - ll^n(0)||,5) ^ 
[B - b,8] which is strictly positive. The result is proved. □ 

Corollary7.3. The sequence of random variables {n~^ Xn) underP^" converges weakly to 
bEIX. 

Proof. Since Sn ^ S and S is strictly increasing and continuous with S(0) - 0, the continuous- 
mapping theorem implies that ^ S~^, see for instance Whitt |18|. Moreover, the 
En's are identically distributed and so £„ ^ E. Since is deterministic we get the 
joint convergence {S~^^,En) ^ (S~^,E) and the continuous-mapping theorem then im- 
plies that S^^iEn) => S~^(i3), since is a continuous function. Since S^^(£'„) = ""^Xn 
and S"'^ (£■) - bEIX this proves the result. □ 
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Corollary7.4. The sequence of random variables [n Xn) underV^ converges weakly to 
E'EI X with E' an exponential random variable independent of E and with parameter a. 

Proof This is a consequence of Corollarv l7.3l together with the fact that n~^x„(0) under 
F]^" converges weakly to E'n by Theorem l6.ll □ 

The above corrolary is simUar to heavy-traffic results for the sojourn time distribu- 
tion ordinary Processor-Sharing queue, see for instance Sengupta |14|, Yashkov |20], 
and Zwart & Boxma [21J. These results may be intuitively explained by the snapshot 
principle, see Reiman [12] : in heavy-traffic conditions the total number of users in the 
system hardly varies over the time scale of a sojourn time. Thus each invididual user 
sees a service rate that is random, determined by the inverse of the total number of 
users in stationarity which has an asymptotically exponential distribution in heavy traf- 
fic, but nearly constant over the duration of its sojourn time. 

It is worth emphasizing that although in the present model the users within each 
of the individual nodes are served in a Processor-Sharing manner, at any given time 
the service rates of users may strongly vary across nodes. Due to the homogenization 
property, however, the empirical distribution of the location of each individual user over 
the course of a long sojourn time in a heavy-traffic regime will be close to the stationary 
distribution n. Hence each individual user will see a :/r-weighted average of the service 
rates in the various nodes, which is only affected by the exponentially distributed total 
number of users in the entire system, just like in an ordinary Processor-Sharing queue. 
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